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(54) Title: DEVICE FOR THE SEGMENTATION OF A DISCRETE SET OF DATA 
(54) Titre: DISPOSITIF DE SEGMENTATION DUN ENSEMBLE DISCRET DE DONNEES 
(57) Abstract 

Device for the segmentation of data corresponding, for example, to image sites, 
for the low rate transmission or storage of said images. The device comprises, firstly, 
an initial segmentation categorizing subunit which issues an initial set of labels. A 
subunit (200) then determines, according to two parallel means (201, 202), conditional 
probabilities, based, in particular, on a propagation operator L, whose parameters are 
defined by given relationships. The probabilities are used in a labelling subunit (300) 
which provides a new set of labels. Finally, a sequencing subunit (400) controls, at 
the input of said subunit (300), the replacement of the initial set of labels with the 
new set and whether the iterative procedure for determination of new labels should 
be continued or discontinued. 

(57) Abreg£ 

Dispositif de segmentation de donn£es correspondant par exemple a des sites 
d' image, pour la transmission et/ou le stockage a faible debit de ces images. Ce 
dispositif comprend tout d'abord un sous-ensemble de segmentation initiale en classes 
deiivrant un jeu de labels initial. Un sous-ensemble (200) opere alors, dans deux 
voies en parallele (201, 202), des determinations de probability conditionnelles, 
notamment a partir d'un opeYateur de propagation L dont les parametres sont d£finis 
par des relations splcifiques. Ces probability sont utilisees dans un sous-ensemble 
d'&iquetage (300) qui ddlivre un nouveau jeu de labels. Enfin, un sous-ensemble 
(400) de sequencement determine a Tentree du sous-ensemble (200) le remplacement 
du jeu de labels initial par le nouveau jeu et la poursuite ou au contraire Tan-fit de ce 
processus iteYatif de determination de nouveaux labels. 
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"DISPOSITIF DE SEGMENTATION D'UN ENSEMBLE DISCRET DE DONNEES" 
Description 

La presente invention concerne un dispositif de 
segmentation d'un ensemble discret de donnees correspondant a 
des sites (m,n) avec m compris entre 1 et M inclus et n entre 1 
et N inclus, a partir de signaux numeriques d f entree 
5 correspondant auxdites donnees. Cette invention trouve des 

applications en. particulier dans le domaine de la segmentation 
d f images, pour realiser par exemple des transmissions de type 
videophonique, ou pour des traitements tels que celui d* images 
fournies par des satellites ou celui d* images obtenues dans le 

10 domaine medical ou en vision robotique, et concerne alors 

toutes sortes d ! images . classiques ou bien a indication de 
distance (entre le capteur de prise de vue et l'objet de la 
scene le plus proche) , et aussi bien des images a distribution 
de niveaux de gris que des images texturees ou des images 

15 associant plusieurs types de structure d' image. 

Tout procede de segmentation d' image a pour objet une 
classification, c'est-a-dire un regroupement des elements 
d f image (ou pixels) en leur affectant une etiquette commune (ou 
label) par region issue de la segmentation. Une telle operation 

20 revient en quelque sorte a construire a partir de donnees 

observables une modelisation facilitant des traitements 
ulterieurs. Cependant, plutot que de construire un modele 
deterministe, on prefere chercher alors a evaluer des 
probabilites attachees aux situations qui sont le reflet de ces 

25 donnees observables, avec une vraisemblance qui est la marque 

du degre de confiance en cette interpretation particuliere . 
Dans 1 1 etablissement d f un tel modele non deterministe, on ne 
peut pas exclure que, dans la classification recherchee, chaque 
pixel soit influence par la classe d ! appartenance de ses 

30 voisins. A partir de cette hypothese, on a songe alors a 

considerer 1' ensemble des pixels d'une image comme un champ de 
Markov ( c 1 est-a-dire comme la generalisation d'une chaine de 
Markov) . 
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La communication "Segmentation of textured images using 
a multiresolution approach " , effectuee par C. Bouman et B. Liu 
lors de la "1988 International Conference on Acoustics, Speech 
and Signal Processing" (tenue a New-York, USA, du 11 au 14 
5 avril 1988) et publiee dans les "Compte-Rendus" de cette 

conference, volume n°2, pp. 1124-1127, presente une methode de 
segmentation d' images qui repose sur une estimation approchee 
du maximum a posteriori d'un champ de Markov aleatoire, a 
partir de I'hypothese qu'une distribution a priori des labels 

10 affectes aux pixels est donnee par une distribution de Gibbs. 

Cette methode precede par des segmentations de plus en plus 
precises, le processus de segmentation avec resolution d 1 image 
croissante etant interrompu lorsque chaque pixel a pu etre 
affecte a une classe. 

15 Le but de 1' invention est de proposer un dispositif de 

segmentation de donnees, et notamment de signaux nuraeriques 
correspondant a des images, qui met en oeuvre une methode de 
segmentation faisant egalement appel a une modelisation par un 
champ de Markov, mais reposant sur une approche differente et 

20 conduisant a un resultat optimal en tres peu d 1 iterations . 

A cet effet, 1* invention concerne un dispositif de 
segmentation plus particulierement caracterise en ce qu'il 
comprend en serie : 

(A) un sous-enserable de segmentation initiale par 

25 extraction d'histogramme, calcul de moyennes sur T classes, et 

classification par recherche de celle des moyennes qui est la 
plus proche du signal calcule ; 

(B) un sous-ensemble de determination et mise a jour de 
probabilites, recevant d'une part lesdits labels et d f autre 

30 part lesdits signaux nuraeriques d 1 entree et comprenant 

lui-raeme : 

(1) en sortie dudit sous-ensemble de segmentation 
initiale, une premiere voie de determination de probabilites 
conditionnelles p(q/a,b,c), comprenant elle-meme en serie : 
35 (a) un circuit de determination de 
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probabilites dites de configuration traduisant l'identite ou la 
non-identite des niveaux correspondant a chacune des quinze 
configurations possibles a l'interieur d 1 un modele 
bidimensionnel en carre (q t a t b,c) de chaine de Markov, q etant 
5 appele dans ledit ensemble a segmenter le label a determiner au 

site dit "enfant" et a f b t c les labels des sites dits 
"parents" ; 

(b) un circuit de calcul de coefficients 
de propagation pour le modele de chaine de Markov ainsi 

10 choisi ; 

(c) un premier circuit de determination de 
la probability conditionnelle p(a/abc) = exp L(q/abc) /Z(abc) , 
ou Z est un facteur de normalisation et ou L(q/abc) est un 
operateur dit de propagation de 1' influence des parents aux 

15 enfants, defini par les expressions (3) a (1*4) : 

L(q/abc) = B(v). d(qa) + B(d) . d(qb) . + B(h). d(qc) 
+ B(v f d). d(qab) + B(d,h) . d(qbc) 
* B(v.h). d(qac) ♦ B(v,d,h). d(qabc) (3) 
exp B(v) = k.S(1123) (4) 
20 exp B(d) = k.S(1213) (5) 

exp B(h) = k.S(123D (6) 
exp B(v,d) = c(v.d) . (exp.B(h) + T-2) .exp (B(v) + B(d)) (7) 
exp B(v,h) * c(v,h) . (exp.B(d) + T-2) .exp (B(v) + B(h)) (8) 
exp B(d.h) = c(d,h).(exp B(v) + T-2) .exp - (B(d) + B(h) ) (9) 



2 c exn B<v d h> = (T-l) , S (1111) /S (1222) 

0 p v exp(B(v) +B(d) +B(h) +B(v,d) +B(v,h) +B(d,h) ) 



(10) 



c(v.d) = S(1112)/(S(1221) ♦ S(1223)) (11) 
c(v,h) = S(1121)/(S(1212) ♦ S(1232)) (12) 
c(d.h) = S(1211)/(S(1122) > S(1233)) (13) 
k = (T-3)/S(123<4) ; (14) 
30 (2) egalement en sortie du sous-ensemble de 

segmentation initiale, et en parallele sur la premiere voie, 
une deuxieme voie de determination de probabilites 
conditionnelles P q (x mn ) t comprenant elle-meme en serie : 

(a) un circuit de calcul de parametres 



BNSDOCID:<WO 9522806A1> 



WO 95/22806 PCT/IB95/00076 



"moyenne" et de parametres "variance" desdites T classes ; 

(b) un deuxieme circuit de determination 
de probabilites conditionnelles ; 

(C) en sortie dudit sous-ensemble de determination et 
5 mise a jour de probabilites, un sous-ensemble d'etiquetage des 

points dudit ensemble de donnees, recevant sur des premiere et 
deuxieme entrees les signaux de sortie desdites premiere et 
deuxieme voies et sur une troisieme entree lesdits signaux 
numeriques d f entree du dispositif, ledit circuit etant prevu 
10 pour delivrer un nouveau jeu de labels qui constitue la mise a 

jour de ladite segmentation initiale ; 

(D) entre les deux sous-ensembles de segmentation 
initiale et de determination et mise a jour de probabilites, un 
sous-ensemble de sequencement, recevant sur des premiere et 

15 deuxieme entrees respectivement le jeu de labels initial obtenu 

au terme de ladite segmentation initiale et le jeu de labels 
resultant de la mise a jour et comprenant en serie : 

(1) un circuit de comparaison desdits jeux de 

labels ; 

20 (2) selon le resultat de ladite comparaison, un 

circuit d* arret du processus de segmentation ou, au contraire, 
pour une nouvelle mise a jour, de substitution dudit jeu de 
labels mis a jour audit jeu de labels initial, ledit jeu de 
labels mis a jour etant envoye vers 1 ! entree desdites premiere 

25 et deuxieme voies. 

Dans un mode particulier de realisation, ledit ensemble 
discret comprend des donnees unimodales et le sous-ensemble de 
segmentation initiale, recevant lesdits signaux d' entree, 
comprend en serie : 

30 (1) un circuit d' extraction de 1 * histogramme des 

niveaux de grandeur exprimes par ces signaux d f entree et qui 
correspondent respectivement a chaque point dudit ensemble ; 

(2) un circuit de calcul , a partir dudit 
histogramme, de moyennes sur T classes correspondant a T 

35 caracteristiques dudit ensemble discret prealablement 
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selectionnees, le nombre predetermine T etant fourni audit 
circuit de calcul ; 

(3) un circuit de classification par recherche, 
pour tous les points dudit ensemble discret, de celle des T 
5 moyennes calculees qui est la plus proche du niveau de chaque 

signal et regroupement correspondant desdits points, ledit 
circuit de classification delivrant des labels qui indiquent 
celle des classes a laquelle appartient chaque point au terme 
de ladite segmentation initiale. 

10 Dans un autre mode de realisation, ledit ensemble 

discret comprend des donnees multimodales pour lesquelles les 
signaux numeriques d* entree correspondant auxdites donnees sont 
regroupes en C canaux, et ledit sous-ensemble de segmentation 
initiale, recevant lesdits signaux d 1 entree, comprend en 

15 serie : 

(1) un etage de segmentation unimodale, recevant 
parmi lesdits signaux d* entree les signaux associes aux M x N 
donnees qui correspondent a 1 1 un des C canaux et comprenant en 
serie : 

20 (a) un circuit d' extraction de 

1 1 histogramme des niveaux de grandeur exprimes par ces signaux 
d* entree et qui correspondent respectivement a chaque point 
dudit ensemble ; 

(b) un circuit de calcul, a partir dudit 
25 histogramme, de moyennes sur T classes correspondant a T 

caracteristiques dudit ensemble discret prealablement 
selectionnees, le nombre predetermine T etant fourni audit 
circuit de calcul ; 

(c) un circuit de classification par 
30 recherche, pour tous les points dudit ensemble discret, de 

celle des T moyennes calculees qui est la plus proche du niveau 
de chaque signal associe a la donnee correspondant a chaque 
point puis par regroupement correspondant desdits points, ledit 
circuit de classification delivrant des labels qui indiquent 
35 celle des classes a laquelle appartient chaque point au terme 
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de la segmentation unimodale ; 

(2) un sous-ensemble de mise a jour de ladite 
segmentation unimodale, comprenant en serie : 

(a) un circuit de mise a jour dudit calcul 
5 de moyennes sur les T classes, recevant d'une part les signaux 

de sortie dudit circuit de classification et d' autre part les 
signaux associes aux M x N x (C-l) autres donnees correspondant 
aux (C-l) autres canaux ; 

(b) un circuit de mise a jour de ladite 
10 classification, delivrant des labels qui indiquent celle des 

classes a laquelle appartient chaque point dudit ensemble 
discret au terme de la segmentation initiale. 

Selon un mode preferentiel de realisation du 
dispositif, les parametres "variance" sont donnes par 
15 l f expression : 

V q.3 = ^.,«<.i . q , (*„.j " M2 «,.j) 2 /Nq 

ou j = 1 a C designe l'indice du canal, x„ , ledit niveau de 

grandeur sur le canal j constate en un site s determine, £(s) 
le label attribue a ce site s, M2 q tj les moyennes sur chaque 

20 canal j independamment , et Nq le nombre total de sites auxquels 

le label q est attribue. 

Selon une variante de realisation du dispositif, 
les parametres "variance" sont les elements de la matrice des 
covariances sur les differents canaux, donnes par 1* expression 

25 suivante : 

V <a,ji,j2 = 2 (s/£(s) « q) ( x s.ji - M2 n.ji)< x s.ji " M2 q . j2 )/ Nc 3 

ou jl,j2 designent les indices des deux canaux concernes, x, . 

S , J 

ledit niveau de grandeur sur le canal j constate en un site s 
determine, £(x) le label attribue a ce site s, M2 q . les 
30 moyennes sur chaque canal j considere independamment, et Nq le 

nombre total de sites auxquels le label q est attribue. 

Dans l'un ou 1' autre cas, une disposition 
avantageuse, dans ce dispositif, consiste a prevoir le renvoi 
de la sortie du sous-ensemble de mise a jour de la segmentation 
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unimodale vers son entree pour realisation d'une nouvelle 
classification, jusqu'a interruption de ce processus iteratif 
selon un critere de convergence determine. 

Quelle que soit la realisation de ce dispositif 
5 qui se trouve concernee, il est en tout cas prevu, dans une 

realisation preferentielle de 1* invention, que le sous-ensemble 
d'etiquetage comprend en serie : 

(1) une memoire des probabilites conditionnelles 
delivrees par lesdites premiere et deuxieme voies et de 

10 probabilites intermediates determinees de fagon interne audit 

sous-ensemble d'etiquetage ; 

(2) un circuit de calcul de grandeurs egales ou 
proportionnelles aux probabilites marginales de chaque label 
pour chaque site, ledit calcul etant conduit selon les etapes 

15 suivantes, pour tous les q t a, b, c possibles : 

(a) on definit F ron (q) comme la probability 

conditionnelle p(q/X mn ) d' avoir au site (m,n) un label q, la 

reference X mn designant 1* ensemble des sites situes au-dessus et 

a gauche du site (m,n) ; 
20 (b) pour chaque b, on examine la valeur 

F = F m _j n _j(b) pour ne la prendre en compte que si elle est 

superieure a un seuil determine ; 

(c) on calcule alors 1* expression : 
H = V B . l . n (a.b).Z 01in . 1 (c.b).p(q/a.b t c).p (| (x IIin )/F 

25 ou : Y m . l n (a f b) = p(a.b/X m . 1( J 

et : Z Btn . 1 (c.b) = Plcb/X,^) 

et l'on effectue une addition cumulee des valeurs de H, pour 
obtenir une grandeur J et une grandeur L m _ x ^(b) ; 

(d) on reprend les etapes (b) et (c) jusqu'a ce 
30 que tous les sites aient ete traites de fagon similaire, la 

valeur cumulee de J etant ajoutee a Y m n (q,c) et a Z m n (q,a) 
puis remise a zero, et la valeur cumulee de L m , 1 n _j(b) 
constituant ladite grandeur egale ou proportionnelle a la 
probability marginale du label b au site (m-l f n-l) ; 
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(e) on reprend les e tapes (b) a (d) jusqu'a ce 
que toutes les combinaisons (q,a,c) aient ete traitees de fagon 
similaire ; 

(3) un circuit d'etiquetage des sites, par comparaison 
5 des valeurs de L m ^ .^(b) et attribution a chaque site (m-l,n-l) 

du label b correspondant a la valeur de L m _ x (b) 

la plus grande. 

Dans le cas de ce dernier mode de realisation, une 

variante avantageuse consiste en ce que, dans le circuit de 
10 calcul de grandeurs egales ou proportionnelles aux probabilites 

marginales de chaque label pour chaque site, il est prevu dans 

l'etape (c) une etape (c 2 ) supplementaire de normalisation par 

division de toutes les valeurs de Y et de 2 par une norme N 

egale, pour tous les b, a 2^ j L (ll . l n _ 1 (b) . 
15 Les particulari tes et avantages de 1' invention 

apparaitront de fagon plus detaillee dans la description qui 

suit et dans les dessins annexes, donnes a titre d 1 exemples non 

limitatifs et dans lesquels : 

- les figures 1 et 4 montrent deux exemples de 
20 realisation du dispositif selon 1' invention ; 

- la figure 2 represente un exemple de chaine de Markov 
d'ordre 3. et la figure 3 montre dans ce cas les differentes 
configurations possibles, dans lesquelles les pixels relies par 
un trait, et seulement ceux-ci, sont identiques. 

25 Avant de decrire plus precisement des exemples 

d' application de 1' invention a la segmentation d ! images, on 
rappelle que, dans une telle application, les images a 
segmenter ont ete prelevees par des moyens qui n'entrent pas 
dans le cadre de cette invention et qui sont directement 

30 dependants de 1 1 application . Les images peuvent provenir de 

satellites, avoir ete obtenues a 1 ' aide d'appareils medicaux, 
etc..-, elles peuvent aussi etre deja le resultat d'un 
pretraitement , comrae par exemple dans le cas d ! images texturees 
ou a chaque image pretraitee correspondent des vecteurs 

35 contenant pour chaque site des caracteristiques de textures 
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tirees prealablement de 1 1 image originale. Dans tous les cas, 
ces images sont disponibles , a 1' entree du dispositif decrit, 
sous forme de signaux nuraeriques. Ces signaux correspondent, 
dans l'exemple decrit en reference a la figure 1, a des donnees 
5 unimodales et, dans l'exemple decrit en reference a la figure 

4, a des donnees multimodales regroupees en C canaux designes 

par j = 1, 2 C. 

Dans la sequence des images successives , chaque image 
numerique est ici definie a 1 ' aide d'un reseau bidimensionnel 
10 de points de coordonnees (x,y) , ou sites s (m,n) avec m = 1 a M 

et n = 1 a N t auquel sont associes des vecteurs dits 
observables I(x,y) quantifies, pour la composante j de chaque 
vecteur, sur V f valeurs distinctes (par exemple de 0 pour le ou 
les points les plus sombres a 255 pour le ou les points les 
15 plus clairs), regulierement espacees en general et appartenant 

a un ensemble C } = g Y , g 2% g 3 g vj . La segmentation de telles 

images consiste a y definir des regions, soit en reduisant le 
nombre des valeurs de quantification des vecteurs observables a 
l'aide d'une technique de transformation qui elimine parmi les 
20 donnees initiales les informations redondantes pour ne 

conserver que les informations essentielles f soit, dans le cas 
de textures, en subdivisant les images en regions spatialement 
coherentes selon certains cri teres d 1 homogeneite . 

Pour respecter le formalisme mathematique qui est en 
25 general adopte lorsqu'on utilise une modelisation par champ de 

Markov (les rappels theoriques concernant cette modelisation 
elle-meme seront effectues plus loin dans la description), on 
dira, dans la suite de la presente description, qu'a chaque 
image X de la sequence des images d f entree est associe un 
30 ensemble S de sites qui correspondent a la localisation de 

chaque pixel. Chaque site i appartient a cet ensemble S qui 
comprend M lignes x N colonnes . 

La segmentation d'une image consiste a associer k 
chaque pixel de celle-ci une etiquette ou label, qu'on notera 
35 2(m,n), m et n designant respec tivement la ligne courante et la 



BNSDOCID:<WO 9522806A1> 



WO 95/22806 




PCT/IB95/00076 



10 

colonne courante du pixel considere, et designant done un site 
determine associe a ce pixel. Ces labels permettent de classer 
chaque pixel dans une region determinee de l f image, sur la base 
de propriete(s) particuliere(s) reconnue(s) comme commune (s) a 
5 la region consideree. II est bien entendu ici que ces labels ne 

sont pas directement observables ou mesurables a partir des 
signaux d' entree representatifs des images, et qu'il faut done 
les determiner. Dans le cadre de la presente invention, le 
nombre de classes, appele par exemple T, doit etre predefini 

10 par 1 'utilisateur , T etant un norabre entier arbitraire. Par 

exemple, dans le cas d'une classification de terrains a partir 
d* images fournies par des satellites, on aura tendance a 
definir T comme etant le nombre de types de terrain differents 
connus a priori. De fa^on similaire, en iraagerie medicale on 

15 choisira le nombre de types de tissu presents dans la partie du 

corps representee par 1' image. 

Dans un premier exemple de realisation, represents sur 
la figure 1 et qui correspond a une segmentation a partir de 
donnees numeriques unimodales ( c 1 es t-a-dire de vecteurs de 

20 donnees observables I(x,y) unimodales), le dispositif de 

segmentation selon 1* invention comprend tout d'abord un sous- 
ensemble de segmentation initiale 100. Ce sous-ensemble 100, 
qui regoit les signaux numeriques representatifs de chaque 
image a segmenter, comprend en serie un circuit 11 d 1 extraction 

25 de 1'histogramme des valeurs observables I(x,y), un circuit 12 

de calcul de moyennes sur T classes correspondant a T 
caracteristiques d f image selectionnees pour la segmentation, et 
un circuit 13 de classification. 

La mise en oeuvre de ce sous-ensemble 100 est la 

30 suivante. En fonction de 1* importance de la segmentation 

recherchee dams 1 ' application consideree, un certain nombre T 
de classes est choisi en correspondance a T caracteristiques 
determinees, par exemple T classes correspondant, dans le cas 
d 1 images en provenance de satellites, a des forets, a des 

35 champs, a des lacs, etc..., ou, dans le cas d ! images medicales. 
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a differents types de tissus corporels (dans la suite de la 
description, on peut prendre par exemple T = 8) . L' histogramme 
des donnees observables l(x,y), extrait a I'aide du circuit 11, 
est fourni au circuit 12 qui le divise en 2T zones de la fagon 
5 suivante : on divise 1 1 histogramme en T zones de surface egale, 
puis on determine dans chacune de ces T zones sa moyenne M1(T). 
Cet ensemble de subdivision revient, en fait, a diviser 
1* histogramme en 2T zones de surface egale deliraitees par 
(2T+1) valeurs observables g(0) , g{l), g(2) g(2T-l), 

10 g(2T), et a prendre les valeurs g(l). g(3) , . . . , g(2T-3) , g(2T- 

1 ) comme valeurs de moyennes . 

Ces T moyennes etant maintenant determinees, le circuit 
13 de classification, qui regoit les signaux de sortie du 
circuit 12 et les (M x N) signaux numeriques d' entree du 

15 dispositif. recherche pour tout pixel (de site i) de 1 1 image la 

moyenne la plus proche. Pour cela, ce circuit 13 effectue la 
determination de la grandeur |g(x) - M1(T)|. Dans cette 
expression g(x) est la valeur quantifiee attachee au pixel 
considere, M1(T) prend successiveraent les T valeurs possibles 

20 (ici, T = 8, chaque pixel etant alors , dans ce cas, range dans 

l'une des T = 8 classes considerees) , et le label correspondant 
a chaque pixel est note ici 2(CL) t ou I'indice CL peut prendre 
T valeurs distinctes. Ce label 2(CL) f determine pour chacun des 
pixels de 1' image, definit une premiere segmentation de 

25 1* image, qu'on appelle ici segmentation initiale. Les M x N 

labels ainsi definis constituent les (M x N) signaux de sortie 
du sous-ensemble de segmentation initiale 100. 

Ce sous-ensemble 100 est suivi d'un sous-ensemble 
200 de raise a jour des probabilites, qui regoit d'une part les 

30 (M x N) signaux d T entree du dispositif et d' autre part les 

(M x N) labels determines par le sous-ensemble 100 et delivre, 
comme on le verra ci-apres,. des probabilites conditionnelles 
p(q/abc) et p q (x) respectivement * Ce sous-ensemble 200 comprend, 
a cet effet, en parallele, des premiere et deuxieme voies 201 

35 et 202 de determination de ces probabilites conditionnelles, et 
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ladite determination est effectuee comrae indique maintenant. 

Si, dans une image, on appelle passe d f un site s 
quelconque 1* ensemble des sites d'indice inferieur (c'est-a- 
dire, par exemple pour un balayage d' image de type television, 
5 1' ensemble constitue par les sites deja balayes sur la meme 

ligne que ce site s et les sites de toutes les lignes 
superieures) , la densite de probability attachee a l'ensemble L 
des labels £ possibles s ! exprime sous la forme d'un produit 
(effectue sur tous les sites) de probabilites d' avoir le label 

10 £(s) au site s = (m,n) , lesdites probabilites etant 

conditionnelles au passe de ce site. La propriete d'un champ de 
Markov dit unilateral se concretise de la fagon suivante : ces 
probabilites conditionnelles (ou probabilites de transition) ne 
dependent pas de tous les labels du passe de ce site, mais 

15 seulement des labels d'un nombre limite des sites appartenant a 

ce passe. Ce nombre de sites limite, qu'on appelle le sous- 
ensemble des parents de s, definit l'ordre du champ de Markov 
ainsi constitue, et le site s lui-raeme est appele par la suite 
le site enfant de ses parents. 

20 Dans les exemples de realisation mentionnes dans la 

presente description, on ne decrira qu f un modele d'ordre 3, 
c'est-a-dire que la probabilite p(2(m,n)) d'obtenir un label 
2(m,n) en un site i(m,n) situe a 1 ' intersection de la m-ieme 
ligne et de la n-ieme colonne ne depend que des probabilites 

25 determinees de fagon similaire pour les sites i(m-l, n), i(m-l, 

n-1), i(m, n-1), avec m et n superieurs a 1. La figure 2 montre 
un tel modele d'ordre 3 : pour simplifier les notations, on a 
appele q le label du site i (ou site "enfant"), et a,b,c les 
labels des trois autres sites "parents". Avec un tel exemple 

30 d'ordre 3. et dans le cas ou toutes les probabilites 

conditionnelles sont presupposees non nulles. 1' expression de 
la probabilite conditionnelle (ou de transition) peut etre 
notee sous forme exponentielle a l*aide d'une fonction (ou 
operateur) L(q/abc). On appelle L un propagateur, car cet 

35 operateur propage une influence des parents aux enfants. On 
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ecrit alors : 

p(q/a.b,c) = exp L(q/a,b.c)/Z(a.b,c) (1) 
Le denominateur Z est simpleraent un facteur de normalisation : 

Z(a,b,c) = ^ exp L(u/a,b,c) , (2) 

u 

5 la sorame 2 sur u etant effectuee sur les T valeurs de labels. 

Ainsi, Z permet d' assurer la con train te de normalisation (on 
doit en effet verifier l'egalite 2 u P(u/a,b,c) = 1, pour tous 
les sites), Dans l'approche ici faite, le propagateur 
L(q/a,b,c) est defini par l'expression (3) suivante : 

10 L(q/a.b.c) = B(v) . d(q f a) + B(d). d(q.b). ♦ B(h) . d(q.c) 

+ B(v.d). d(q.a.b) ♦ B(d f h). d(q.b.c) 
+ B(v,h). d(q.a.c) + B(v.d.h). d(q.a.b.c) (3) 
Dans cette expression, v, d, h designent respectivement le site 
parent vertical, diagonal ou horizontal du site enfant i, les 

15 B(.) sont des pararaetres independants constituant des 

coefficients dits de propagation, et les d(.) sont egaux a 1 si 
les labels correspondant aux indices entre parentheses sont 
identiques, ou a 0 sinon. Dans le propagateur ainsi defini, les 
sept valeurs B(.), qui representent les parametres du modele 

20 markovien, caracterisent la tendance d'un site enfant de se 

voir decerner le meme label que ses parents : done, si par 
exeraple la constante B(v) est positive et grande , le site 
enfant aura forte tendance a avoir le meme label que son parent 
vertical . 

25 La connaissance de la probability p(q/a,b,c) depend 

done de celle du propagateur L, e'est-a-dire de celle des 
parametres B(.)* On peut alors demontrer que la connaissance de 
ces parametres depend elle-meme de celle des probabilites de 
configuration a l'interieur du modele de chaine de Markov 

30 choisi (celui represents sur la figure 2). Ces probabilites de 

configuration sont no tees respectivement S ( 123*0 t S(1233), 
S(1232). S(123D. S(1223), S(1222) f S(1221) . S (1213). S(1212). 
S(1211), S(1123). S(1122), S(1121), S(1112). S(llll). La 
signification de ces grandeurs S(.) dans le cas du modele de la 
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figure 2 est visualisee sur la figure 3» qui associe a chaque 
probability de configuration S(.) la representation de la 
configuration correspondante f a savoir celle pour laquelle les 
pixels relies par un trait, et seulement ceux-ci, sont 
5 identiques. Ces probabilites de configuration sont determinees 

par simple comptage du nombre de fois ou, dans 1' image 
concernee, apparait une configuration determinee. 

Ce comptage est realise dans la premiere voie 201 qui 
comprend a cet effet un circuit 211 de determination de ces 

10 probabilites de configuration S(.) dans 1' image segmentee 

presente en sortie du sous-ensemble 510. Ce circuit 211 est 
suivi d'un circuit 212 de calcul des coefficients de 
propagation B(.)- On a vu plus haut que la connaissance de ces 
coefficients dependait de celle des probabilites S(.). Sans 

15 expliciter l f ensemble des raisonnements qui permettent 

d'etablir cette dependance, on ne donnera ici que les relations 
(4) a (lU) qui lient les B(.) et les S(.), et qui, dans le cas 
ou le nombre T est superieur ou egal a 4, sont les suivantes : 
exp B(v) = k.S(1123) (4) 

20 exp B(d) = k.S(1213) (5) 

exp B(h) = k.S(123D (6) 

exp B(v,d) = c(v.d) . (exp.B(h) * T-2).exp (B(v) ♦ B(d)) (7) 

exp B(v B h) = c(v,h) . (exp.B(d) + T-2) .exp (B(v) + B(h)) (8) 

exp B(d,h) = c(d,h).<exp B{v) ♦ T-2). exp (B(d) + B(h)) (9) 

25 exp B(v d h)= (T-l) .S(llll)/S(1222) ( , 

° P ' ' exp(B(v) +B(d) +B<h) +B(v,d) +B(v,h) +B(d,h) ) 

c(v.d) = S(1112)/(S(1221) ♦ S(1223)) (11) 
c(v.h) = S(1121)/(S(1212) + S(1232)) (12) 
c(d.h) = S(1211)/(S(1122) + S(1233)> (13) 
k = (T-3)/S(123^). (14) 
30 Dans le cas ou T = 3 . les equations (4) a (14) res tent 

verifiees, mais la valeur de k est arbitraire (toute valeur de 
k donne lieu a la meme valeur de L(q/a,b,c) et conduit done a 
des segmentations identiques) . On peut par exemple choisir k 
tel que B(v) + B(d) + B(h) = 0, ce qui implique que k soit egal 
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a la racine cubique de 1/S(1123) x S(1213) x S(1231), mais tout 
autre choix de k convient aussi . De meme, dans le cas ou 
T = 2, toute valeur de k donne lieu aux memes valeurs 
L(q/a,b,c), et done a des segmentations identiques. Les 
5 coefficients B(.) etant ainsi determines, un circuit 213 prevu 

en serie avec les circuits 211 et 212 determine, conformement a 
l'expression (1), la probability conditionnelle (ou de 
transition) p(q/a,b,c) . 

La deuxieme voie 202 regoit, comme la premiere, les 

10 (M x N) labels de sortie du sous-ensemble 100, et les fournit a 

un circuit 221 de calcul de parametres de classe (ici, moyennes 
et variances). En effet, le circuit 12 avait effectue un 
premier calcul de moyennes sur T classes, mais ces classes 
avaient ete, pour 1 1 initialisation du f onctionnement du 

15 dispositif de segmentation, constitutes de fagon arbitraire, 

par division de 1 1 histogramme des valeurs observables en T 
zones. Cette segmentation initiale arbitraire etant maintenant 
remplacee en sortie du sous-ensemble 100 par une segmentation 
repondant desormais au moins approximativement a un critere de 

20 classement, il est done possible d'effectuer un nouveau calcul 

de parametres pour chacune des T classes repondant aux T labels 
possibles pour les pixels, plus precisement un nouveau calcul 
de la moyenne, notee M2(T). ainsi que le calcul de la variance 
(elevee au carre) correspondante , notee V2(T) . Les T valeurs 

25 M2(T) et V2(T) ainsi calculees sont les 2T signaux de sortie du 

circuit 221. 

Ces 2T signaux sont fournis a un circuit 222 de 
determination de la probability conditionnelle d'un niveau de 
gris x pour un label q constate en un site i. On sait en effet 
30 que la segmentation finale doit etre en accord avec les donnees 

initiales, qui resultent d'une observation ou d'une mesure. 
Dans le cas de donnees observables I(x,y) unimodales, cette 
probability conditionnelle, notee p q (x), est donnee pour chaque 
classe par 1* expression (15) suivante : 
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p q (x) - i— -exp <x-M2<T))» (15) 

q V2(T)v/2tc 2V2(T) 2 

Les probabilites conditionnelles de sortie des deux 
voies 201 et 202, presentes en sortie des circuits 213 et 222, 
constituent les deux types de signaux de sortie du sous- 
5 ensemble 200 de mise a jour des probabilites, et sont alors 

envoyees vers deux entrees correspondantes d'un sous-ensemble 
300 d'etiquetage des pixels. Ce sous-ensemble 300 regoit 
egalement, sur une troisierae entree, les (M x N) signaux 
numeriques d f entree du dispositif, et determine des 

10 probabilites dites marginales, a partir desquelles sont a leur 

tour determines de nouveaux labels £(.) definissant une 
nouvelle segmentation de 1' image. Ces probabilites marginales 
sont determinees comme indique ci-dessous, dans la description 
du f onctionnement du sous-ensemble 300. 

15 On rappellera cependant, prealablement , quelques 

notations. Le site qui est momentanement traite est le site 
(ra,n) , avec m compris entre 1 et M inclus et n entre 1 et N 
inclus. La partie d' image qui est situee au-dessus et a gauche 
du site (m,n) est notee X nin . La probability conditionnelle 

20 P(q/X mn ) d'avoir au site (m t n) le label q est notee F mn (q) . On 

designe enfin par p(q,a,b,c/X n)n ) la probability marginale 
recherchee d f avoir les labels (ou etiquettes) 2„ = q f £ . = a 

^m-i.n-i ~ b » c ni. n-i = c sur les sites correspondants . 

La determination de p(q,a,b,c/X mn ) pour chaque site est 

25 realisee, dans le sous-ensemble 300. comme indique maintenant, 

les donnees utilisees par ce sous-ensemble etant les 
suivantes : les probabilites conditionnelles p q (x) delivrees par 
la voie 202 (pour lesquelles 1' ensemble des x pour chaque site 
m,n correspond a I 1 ensemble des donnees observables I), les 

30 probabilites conditionnelles p(q/a,b,c) delivrees par la voie 

201, ainsi que les probabilites suivantes : 
Y-».„<a.b) = P(a.b/X.. liM ) et ^.....(cb). = plc.b/X.,^), 
determinees pour tous les q, a, b, c possibles et de fagon 
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interne audit sous-ensemble 300. Toutes ces probabilities sont 
stockees dans une memoire 301. 

Un circuit de calcul 302 situe en sortie de cette 
memoire precede alors, pour tous les q. a, b, c possibles, aux 
5 determinations suivantes, par etapes successivement enumerees : 

(1) pour chaque b t on examine la valeur 

F " F m-i. n -i( b ) • n'est pas prise en compte si elle est 

inferieure a un seuil determine, par exemple 1CT 8 ; 

(2) si au contraire cette valeur depasse ledit seuil, 

10 on calcule une grandeur H, notee H nm (q , a, b , c) selon 1' expression 

(16) : 

H = Y m _i.nU.b).Z ni . n . 1 (c,b).p(q/a.b ( c).p n (x mn )/F (16) 
F venant d'etre definie, et les Y et Z ayant ete prealableraent 
initialises a zero (puis calcules de fagon particuliere 
15 -explicitee plus loin- pour la premiere ligne et pour la 

premiere colonne) ; 

(3) cette valeur de H est additionnee a une grandeur J, 
initialisee a zero (cette addition est en fait une addition 
cumulee des valeurs de H) et a une grandeur L,,,^ ..^(b) 

20 (egalement initialisee a zero) egale ou proportionnelle a la 

probability marginale du label b au site (m-l.n-1) et donnee 
par 1* expression (17) : 

L m -i.n.i(b) = 2 qnc H nm (q.a.b.c) (17) 
puis on reprend I'etape (1) tant que tous les b n'ont pas ete 

25 traites de fagon similaire ; 

(4) la valeur ainsi obtenue pour J est ajoutee a 
Y m .„(q.c) et a Z m>(l (q,a) ; 

(5) apres remise de J a zero, on reprend les etapes (1) 
a (4) tant que toutes les combinaisons (q,a,c) n'ont pas ete 

30 traitees de fagon similaire. 

Lorsque toutes les combinaisons (q.a.c) ont ete 
ef fectivement prises en compte, le circuit de calcul 302 
delivre pour toutes ces combinaisons les grandeurs L m-1 nM (b) t 
Y m . n (qc), Z mn (qa). Pour eviter que les variables conservees en 
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memoire ne deviennent trop petites pour un traitement numerique 
precis, on procede a une operation de normalisation, qui 

consiste dans le cas present a determiner une norme N definie 

comme la somme des grandeurs L^j ^tb) pour tous les b, et a 
5 diviser toutes les grandeurs Y m n (qc) et 2 m n (qa) par cette norme 

N ra n . On renvoie ensuite du circuit 302 vers la memoire 301 les 
valeurs normalisees qui sont necessaires pour les calculs 
ulterieurs. Comme on l ! a vu plus haut, les probabilites Y et Z 
sont utilisees de fagon interne au sous-ensemble 300 , pour le 
10 calcul de H selon l'expression (16), tandis que les valeurs de 

^m-i.n-i (b) constituent les signaux de sortie du circuit de 
calcul 302. 

Ce circuit 302, qui a done calcule des grandeurs 
egales ou proportionnelles aux probabilites marginales de 

15 chaque label pour chaque site, est alors suivi d' un circuit 303 

d' etiquetage des sites (sauf les sites des premiere et derniere 
lignes et des premiere et derniere colonnes , qui sont traites 
de maniere particuliere comme on le decrit ci-dessous) . Cet 
etiquetage des sites est realise par comparaison des valeurs de 

20 L^.j n .j(b) et attribution au site (m-1, n-1) du label b 

correspondant a la valeur de L la plus grande. 

Dans le cas d'un site de la premiere ligne ou de la 
premiere colonne, on n'effectue pas d ' etiquetage , mais 
simplement on reraplit progressivement la memoire 301 (comme 

25 indique ci-apres) avec les valeurs necessaires pour le 

traitement des lignes et colonnes suivantes. Pour le premier 
site (1,1), la probability de distribution des labels est 
donnee, apres normalisation, par : 

F K1 (q) = P^,.,)/^^,,) (18) 

30 expression dans laquelle p^Xj J est egale pour tous les 

labels. Cette valeur F x #1 (q) est stockee dans la memoire 301 . 

Pour le traitement des sites de la premiere ligne, on 
calcule les F Kn (q) a partir des expressions (19) a (22) 
suivantes : 
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F' l( „(q) = S c Y l#n {q f c) (19) 
ou : F' I>n (q) = S c F^^lc). P hor (q/c) . P q (x l-n ) (20) 
expression (20) dans laquelle P hor (q/c) f calcule dans la 
premiere voie 201, est donne par : 
5 P hor (q/c) = exp((B(h).d(q t c))-log(exp.B(h)*T-l)) (21) 

Ce calcul est, comrae precedemment , suivi d f une normalisation 
(division par N 1 n = 2 q F\ n (q)), qui conduit finalement a : 
F l<n (q) = F\. n (q>/^ F\. n (q) (22) 
Pour le traiteraent des sites de la premiere colonne, on 
10 calcule de fagon similaire les F m j(q) a partir des expressions 

(23) a (26) suivantes : 

FY^q) = S, Z BBl (q f a) (23) 
ou : F' mtl (q) = 2, F niM . 1 ( a) . P^ r (q/a) . Pi) ( x lu . l ) (24) 
expression (24) dans laquelle P ver (q/a) , egalement calcule dans 
15 la premiere voie 201, est donne par : 

P ver (q/a) = exp((B(v).d(q.a)) - log(exp.B(v) +T-1) ) (25) 
Ce calcul est, la encore, suivi d'une normalisation conduisant 
finalement a : 

F-.i(q) = F'-.itqJA, F'-.i(q) (26) 

20 Les valeurs de Y, Z t F sont au fur et a mesure conservees dans 

la memoire 301, pour permettre les calculs ulterieurs. 

Dans le cas d'un site de la derniere colonne (n = N) , 
le circuit de calcul 302 evalue L Bl _ 1<n (a) = 5^ z m n (qa) et le 
label qui rend L m _ Kll maximal est attribue au site (m-1, N) . De 

25 merae, dans le cas d'un site de la derniere ligne (ra = M) f le 

circuit de calcul 302 evalue L^.^c) = 2^ Y m n (qc) et le label 
qui rend L m n .j maximal est attribue au site (M f n-1). Pour le 
site (M,N), le circuit 302 evalue L m n (q) = F m . n (q) et le label 
qui rend L m n maximal est attribue au site (M,N). L f ensemble des 

30 labels ainsi delivres par le circuit 303 forme un jeu de labels 

modifie constituant une mise a jour du jeu de labels initial. 

Un sous-ensemble de sequencement 400 re^oit alors, sur 
des premiere et deuxierae entrees, ledit jeu de labels initial 
et le jeu de labels modifie resultant de la mise a jour, et en 
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effectue tout d'abord la comparaison dans un circuit de 
comparaison 410. Selon que les jeux de labels different ou non 
pour plus d'un nombre determine de pixels (par exemple 40) f un 
circuit 420 situe a la suite du circuit 1*10 commande 1' arret du 
5 processus de segmentation (et le jeu de labels mis a jour est 

definitivement selectionne) ou au contraire la substitution, a 
1' entree des deux voies 201 et 202, dudit jeu de labels mis a 
jour au jeu de labels initial. S'il y a substitution, une 
nouvelle raise en oeuvre des sous-ensembles 200 a kOO 
10 intervient. S'il y a au contraire arret de ce processus, la 

mise en oeuvre du dispositif de segmentation selon l 1 invention 
est terminee, et le jeu de labels ainsi selectionne constitue 
les signaux de sortie definitifs de ce dispositif de 
segmentation . 

15 Bien entendu, 1' invention n'est pas limitee a cet 

exemple de realisation, a partir duquel des variantes peuvent 
etre proposees sans sortir pour cela du cadre de cette 
invention. En particulier, dans un deuxieme exemple de 
realisation, represents sur la figure ^ et qui correspond a une 

20 segmentation a partir de donnees nuraeriques multiraodales 

(c 'est-a-dire de vecteurs de donnees observables I(x,y) 
multimodales ) , le dispositif de segmentation selon 1* invention 
comprend cette fois tout d'abord un sous-ensemble de 
segmentation initiale 510- Ce sous-ensemble 510, qui regoit 

25 lesdites donnees numeriques, comprend un etage 101 de 

segmentation unimodale recevant parmi les M x N x C donnees 
d* entree les M x N donnees qui correspondent a un seul des 
canaux. Cet etage 101 comprend lui-meme en serie exactement les 
memes elements que le sous-ensemble de segmentation initiale 

30 100, a savoir un circuit 11 d 1 extraction de 1 1 histogramme 

desdites M x N donnees, un circuit 12 de calcul de moyennes sur 
T classes correspondant a T caracteristiques d' image 
selectionnees pour la segmentation, et un circuit 13 de 
classification ♦ La mise en oeuvre de cet etage 101 est la meme 

35 que celle du sous-ensemble 100 de la figure 1 et les M x N 
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labels £(CL) determines, qui constituent les (M x N) signaux de 
sortie de l'etage de segmentation 101. definissent une premiere 
segmentation unimodale de 1* image, 

Dans une premiere realisation, cette segmentation 
constitue la sortie de l'etage 510 (cette situation est 
representee sur la figure 4 a 1 ' aide de la connexion en trait 
discontinu en sortie du circuit 13). Dans une deuxieme 
realisation, l'etage 101 est suivi en serie d'un sous-ensemble 
de mise a jour de cette segmentation unimodale, comprenant en 
serie un circuit 51^ de mise a jour du calcul de moyennes 
effectue par le circuit 12, et un circuit 515 de mise a jour de 
la classification effectuee par le circuit 13 . Le circuit 5l4 # 
identique au circuit 221 decrit plus loin, regoit d'une part 
les signaux de sortie du circuit de classification 13 et 
d* autre part les signaux associes aux M x N x C donnees 
d' entree. Le circuit 515 qui le suit delivre des labels qui 
indiquent celle des classes a laquelle appartient chaque point 
au terme de ladite segmentation initiale. Dans ladite deuxieme 
realisation, cette segmentation constitue la sortie de l'etage 
510. 

Une troxsieme realisation de l'etage 510, fonctionnant 
cette fois de fagon iterative, prevoit, comme indique sur la 
figure 4, d'envoyer la sortie du sous-ensemble de mise a jour 
de la segmentation unimodale vers son entree ( c ' es t-a-dire la 
sortie du circuit 515 vers 1' entree du circuit 514). Ainsi, 
apres calcul des valeurs moyennes de chaque classe separement 
(par la prise en compte des pixels appartenant a cette classe, 
a l'aide de formules equivalentes a 1* expression (15)), puis 
classement des pixels de 1' image suivant les C canaux par un 
critere tel que celui d'une distance minimale, il est possible 
d'obtenir une nouvelle classification des pixels. Ce processus 
peut etre repete, jusqu*a convergence vers un minimum de la 
distance par exemple. Cette segmentation constitue la sortie de 
l'etage 510 dans sa troisieme realisation. 

Le sous-ensemble de segmentation initiale 510 est 
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suivi, en serie, du sous-ensemble 200 de mise a jour de 
probabilites , qui regoit comme precedemraent les signaux 
d 1 entree du dispositif et les (M x N) labels determines par le 
sous-ensemble 510 et delivre les probabilites conditionnelles 
5 p(q/abc) et p q (x) . Le sous-ensemble 200 comprend, comme dans le 

cas de la figure 1, les premiere et deuxieme voies 201 et 202, 
identiques a la realisation deja decrite, sauf en ce qui 
concerne le circuit 221 de calcul de parametres de classe. 

Dans une premiere variante de realisation, le 

10 circuit 221 calcule en effet les moyennes et les variances, sur 

chaque canal independamment . Les moyennes, notees M2 q j# sont 
donnees par 1' expression (27) suivantes : 

M2 m.J = 2 (*/«<*> . n > x S!j /Nq (27) 
ou j = 1 a C designe I'indice du canal, x s#j le niveau de 

15 grandeur sur le canal j constate en un site s determine, fi(s) 

le label attribue a ce site s, et Nq le nombre total de sites 
auquel le label q est attribue. Les variances V sont, elles, 
donnees ici par l 1 expression (28) suivante : 
v q 2 j = Z fs/fi(s) s q) (x s j - M2 q<j ) 2 /Nq (28) 

20 Dans une autre variante de realisation, le jeu de parametres de 

classe comprend, en plus des moyennes M2 q j( les elements de la 
matrice des covariances sur les differents canaux, donnes par 
1* expression (29) suivante : 

V q,ji.j2 - 2 (s/ C (s) . M)( x s.jrM2 q(jl )(x s>jl -M2 q j2 )/Nq (29) 
25 Dans l ! une ou 1* autre de ces variantes, les valeurs M2 et V 

ainsi calculees sont, comme precedemment t les 2T signaux de 

sortie du circuit 221. 

Ces signaux sont fournis au circuit 222 de 

determination de probability conditionnelle . Dans le cas de la 
30 premiere variante de realisation du circuit 221, cette 

probability conditionnelle p q (x) est donnee pour chaque classe 

par 1' expression (30) suivante : 

P q (x> - ^_ —3. D exp (-0,5 j a g ' j ) (30) 
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Dans le cas de la deuxieme variante de realisation du circuit 
221, la probability conditionnelle determinee par le circuit 
222 est cette fois donnee, pour chaque classe, par 1' expression 
(31) suivante : 

1 J ' D 



expression (31) dans laquelle b q jlj2 est la ma trice inverse de 
V q jlj2 . Dans l'un et I 1 autre cas, les probabilites 
conditionnelles de sortie des deux voies 201 et 202, presentes 
en sortie des circuits 213 et 222, constituent les signaux de 
10 sortie du sous-ensemble 200 envoyes vers les deux entrees du 

sous-ensemble 300 d'etiquetage des pixels. Ce sous-ensemble 300 
et le sous-ensemble de sequencement ^00 ne sont pas modifies 
par rapport a la realisation de la figure 1. 
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REVENDICATIONS : 

1. Dispositif de segmentation d'un ensemble discret de 

donnees correspondant a des sites (m,n) avec m compris entre 1 
et M inclus et n entre 1 et N inclus, a partir de signaux 
5 numeriques d* entree correspondant auxdites donnees, caracterise 

en ce qu'il comprend en serie : 

(A) un sous-ensemble de segmentation initiale par 
extraction d ' histogramme , calcul de raoyennes sur T classes, et 
classification par recherche de celle des moyennes qui est la 

10 plus proche du signal calcule ; 

(B) un sous-ensemble de determination et raise a jour de 
probabilites, recevant d'une part lesdits labels et d ! autre 
part lesdits signaux numeriques d 1 entree et coraprenant 
lui-raeme : 

15 (1) en sortie dudit sous-ensemble de segmentation 

initiale, une premiere voie de determination de probabilites 
conditionnelles p(q/a,b,c), comprenant elle-meme en serie : 

(a) un circuit de determination de 
probabilites dites de configuration traduisant l'identite ou la 

20 non-identite des niveaux correspondant a chacune des quinze 

configurations possibles a l'interieur d'un modele 
bidimensionnel en carre (q,a,b.c) de chaine de Markov, q etant 
appele dans ledit ensemble a segmenter le label a determiner au 
site dit "enfant" et a,b,c les labels des sites dits 

25 "parents" ; 

(b) un circuit de calcul de coefficients de 
propagation pour le modele de chaine de Markov ainsi choisi ; 

(c) un premier circuit de determination de la 
probability conditionnelle p(a/abc) = exp L(q/abc) /Z(abc) , ou Z 

30 est un facteur de normalisation et ou L(q/abc) est un operateur 

dit de propagation de 1* influence des parents aux enfants, 
defini par les expressions (3) a (14) : 

L(q/abc) = B(v). d(qa) ♦ B(d). d(qb) . + B(h) . d(qc) 
♦ B{v,d). d(qab) * B(d.h) . d(qbc) 
35 + B(v.h). d(qac) * B(v.d.h). d(qabc) (3) 
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exp B(v) = k.S(1123) (A) 

exp B(d) = k.S(1213) (5) 

exp B(h) = k.S(123D (6) 

exp B(v,d) = c(v,d) . (exp.B(h) + T-2).exp (B(v) + B(d)) (7) 

exp B(v t h) = c(v,h) . (exp.B(d) + T-2).exp (B(v) + B(h)) (8) 

exp B(d.h) = c(d,h).(exp B(v) + T-2).exp (B(d) * B(h)) (9) 



extl r/ v d hi= (T-l) .S<llll)/S(1222) 

* y ' U/ exp(B(v) +B(d) +B(h) +B(v,d> +B(v,h) +B(d,h> ) 



(10) 



c(v.d) = S(1112)/(S(1221) ♦ S(1223)) (11) 
c(v.h) = S(1121)/(S(1212) ♦ S(1232)) (12) 
10 c(d.h) = S(1211)/(S(1122) + S(1233)) (13) 

k = (T-3)/S{i23*») ; (14) 

(2) egalement en sortie du sous-ensemble de 
segmentation initiale, et en parallele sur la premiere voie, 
une deuxieme voie de determination de probabilites 
15 conditionnelles P q (x mn ), comprenant elle-meme en serie : 

(a) un circuit de calcul de pararaetres 
"raoyenne" et de parametres "variance" desdites T classes ; 

(b) un deuxieme circuit de determination de 
probabilites conditionnelles ; 

20 (C) en sortie dudit sous-ensemble de determination et 

mise a jour de probabilites, un sous-ensemble d f etiquetage des 
points dudit ensemble de donnees, recevant sur des premiere et 
deuxieme entrees les signaux de sortie desdites premiere et 
deuxieme voies et sur une troisieme entree lesdits signaux 

25 numeriques d' entree du dispositif, ledit circuit etant prevu 

pour delivrer un nouveau jeu de labels qui constitue la mise a 
jour de ladite segmentation initiale ; 

(D) entre les deux sous-ensembles de segmentation 
initiale et de determination et mise a jour de probabilites, un 

30 sous-ensemble de sequencement , recevant sur des premiere et 

deuxieme entrees respectivement le jeu de labels initial obtenu 
au terme de ladite segmentation initiale et le jeu de labels 
resultant de la mise a jour et comprenant en serie : 

(1) un circuit de comparaison desdits jeux de 
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labels ; 

(2) selon le resultat de ladite comparaison, un 
circuit d f arret du processus de segmentation ou, au contraire, 
pour une nouvelle mise a jour, de substitution dudit jeu de 
5 labels mis a jour audit jeu de labels initial, ledit jeu de 

labels mis a jour etant envoye vers 1* entree desdites premiere 
et deuxieme voies. 

2. Dispositif de segmentation selon la revendication 
1, caracterise en ce que ledit ensemble discret comprend des 

10 donnees unimodales et en ce que le sous-ensemble de 

segmentation initiale, recevant lesdits signaux d ! entree, 
comprend en serie : 

(1) un circuit d' extraction de 1 1 his togramme des 
niveaux de grandeur exprimes par ces signaux d' entree et qui 

15 correspondent respectivement a chaque point dudit 

ensemble ; 

(2) un circuit de calcul, a partir dudit 
histogramme, de moyennes sur T classes correspondant a T 
caracteristiques dudit ensemble discret prealablement 

20 selectionnees , le nombre predetermine T etant fourni audit 

circuit de calcul ; 

(3) un circuit de classification par recherche, 
pour tous les points dudit ensemble discret, de celle des T 
moyennes calculees qui est la plus proche du niveau de chaque 

25 signal et regroupement correspondant desdits points, ledit 

circuit de classification delivrant des labels qui indiquent 
celle des classes a laquelle appartient chaque point au terme 
de ladite segmentation initiale. 

3. Dispositif de segmentation selon la revendication 
30 1, caracterise en ce que ledit ensemble discret comprend des 

donnees multimodales pour lesquelles les signaux nuraeriques 
d f entree correspondant auxdites donnees sont regroupes en C 
canaux et en ce que ledit sous-ensemble de segmentation 
initiale, recevant lesdits signaux d' entree, comprend en 
35 serie : 
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(1) un etage de segmentation unimodale, recevant 
parmi lesdits signaux d' entree les signaux associes aux M x N 
donnees qui correspondent a l'un des C canaux et comprenant en 
serie : 

5 (a) un circuit d ! extraction de l ! histogramme 

des niveaux de grandeur exprimes par ces signaux d f entree et 
qui correspondent respectivement a chaque point dudit 
ensemble ; 

(b) un circuit de calcul, a partir dudit 
10 histogramme , de raoyennes sur T classes correspondant a T 

caracteristiques dudit ensemble discret prealablement 
selectionnees , le nombre predetermine T etant fourni audit 
circuit de calcul ; 

(c) un circuit de classification par 

15 recherche, pour tous les points dudit ensemble discret, de 

celle des T moyennes calculees qui est la plus proche du niveau 
de chaque signal associe a la donnee correspondant a chaque 
point puis par regroupement correspondant desdits points, ledit 
circuit de classification delivrant des labels qui indiquent 

20 celle des classes a laquelle appartient chaque point au terme 

de la segmentation unimodale ; 

(2) un sous-ensemble de mise a jour de ladite 
segmentation unimodale, comprenant en serie : 

(a) un circuit de mise a jour dudit calcul de 
25 moyennes sur les T classes, recevant d f une part les signaux de 

sortie dudit circuit de classification et d' autre part les 
signaux associes aux M x N x (C-l) autres donnees correspondant 
aux (C-l) autres canaux ; 

(b) un circuit de mise a jour de ladite 

30 classification, delivrant des labels qui indiquent celle des 

classes a laquelle appartient chaque point dudit ensemble 
discret au terme de la segmentation initiale. 

4. Dispositif de segmentation selon la revendication 

3, caracterise en ce que les parametres "variance" sont donnes 
35 par 1' expression : 
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V Q.i 3 2 (s/ fitS ) = .,) <*..j - M2 qtj ) 2 /Nq 

ou j = 1 a C designe I'indice du canal, x s tj ledit niveau de 
grandeur sur le canal j constate en un site s determine, £(s) 
le label attribue a ce site s, M2 q . les moyennes sur chaque 
5 canal j independamment , et Nq le nombre total de sites auxquels 

le label q est attribue. 

5- Dispositif de segmentation selon la revendication 

3, caracterise en ce que les parametres "variance" sont les 
elements de la matrice des covariances sur les differents 

10 canaux, donnes par 1' expression suivante : 

V q,ji.j2 = S (s/*(s) . q) (*s.ji - M 2 qijl )(x K<jl - M2 q j2 )/Nq 
ou jl,j2 designent les indices des deux canaux concernes , x t . 
ledit niveau de grandeur sur le canal j constate en un site s 
determine, £(x) le label attribue a ce site s, M2 q s les 

15 moyennes sur chaque canal j considere independamment, et Nq le 

nombre total de sites auxquels le label q est attribue. 
6. Dispositif de segmentation selon l f une des 

revendications 4 et 5. caracterise en ce que la sortie du sous- 
ensemble de mise a jour de la segmentation unimodale est 

20 renvoyee vers son entree pour realisation d'une nouvelle 

classification, jusqu'a interruption de ce processus iteratif 
selon un critere de convergence determine. 
7* Dispositif de segmentation selon l'une des 

revendications 1 a 6, caracterise en ce que le sous-ensemble 

25 d'etiquetage comprend en serie : 

(1) une meraoire des probabilites conditionnelles 
delivrees par lesdites premiere et deuxieme voies et de 
probabilites intermediaires determinees de fagon interne audit 
sous-ensemble d'etiquetage ; 

30 (2) un circuit de calcul de grandeurs egales ou 

proper tionnelles aux probabilites marginales de chaque label 
pour chaque site, ledit calcul etant conduit selon les etapes 
suivantes, pour tous les q, a, b, c possibles : 

(a) on definit F ni n (q) comme la probability 
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conditionnelle p(q/X imi ) d'avoir au site (m t n) un label q, la 
reference X mn designant l 1 ensemble des sites situes au-dessus et 
a gauche du site (m,n) ; 

(b) pour chaque b f on examine la valeur 
5 F = F m-i. n -i( b ) pour ne la prendre en compte que si elle est 

superieure a un seuil determine ; 

(c) on calcule alors l 1 expression : 
H « Y m -i.n(a.b).Z m . n . 1 (c f b).p(q/a.b # c).p q (x mn )/F 

ou : Y m . Kn (a.b) = p{a.b/X B . Kn ) 
10 et : Z llltH (c,b) = P(c.b/X 1)ltn . 1 ) 

et l f on effectue une addition cumulee des valeurs de H, pour 
obtenir une grandeur J et une grandeur L m-1 n _j(b) ; 

(d) on reprend les etapes (b) et (c) jusqu'a 
ce que tous les sites aient ete traites de fagon similaire, la 

15 valeur cumulee de J etant ajoutee a Y m n (q,c) et a Z m n (q t a) puis 

remise a zero, et la valeur cumulee de 

L m _ 1 .n-xfb) constituent ladite grandeur egale ou proportionnelle 
a la probabilite marginale du label b au site (m-l,n-l) ; 

(e) on reprend les etapes (b) a (d) jusqu'a 
20 ce que toutes les combinaisons (q,a,c) aient ete traitees de 

faQon similaire ; 

(3) un circuit d'etiquetage des sites, par comparaison 
des valeurs de L ni _ L n-1 (b) et attribution a chaque site (m-l ( n-l) 
du label b correspondant a la valeur de 

25 L m-i.ii-i< b ) la P lus S^ande. 

8. Dispositif de segmentation selon la revendication 

7. caracterise en ce que, dans le circuit de calcul de 
grandeurs egales ou proportionnelles aux probabilites 
marginales de chaque label pour chaque site, il est prevu dans 

30 l'etape (c) une etape (c 2 ) suppleraentaire de normalisation par 

division de toutes les valeurs de Y et de Z par une norme N m n 
egale, pour tous les b, a 2^L m _ l<n ., (b) . 
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